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ABSTRACT 



Aims. In this paper we study the influence of the ionizing cluster mass on the emission line spectrum of Hn regions in order to 
determine the influence of low mass clusters on the integrated emission line spectra of galaxies. 

Methods. For this purpose, we present a grid of photoionization models that covers metallicities from Z = 0.001 to Z = 0.040, ages 
from 0. 1 to 10 Ma (with a time step of 0. 1 Ma), and cluster initial masses from 1 to 10 7 M Q . The stellar masses follow a Salpeter initial 
mass function (IMF) in an instantaneous burst mode of star formation. We obtain power-law scale-relations between emission-line 
luminosities and ionizing cluster masses from the grids and we evaluate the dependences on the ionizing cluster mass for some line 
luminosities, equivalent widths and line ratios. 

Results. Power-law scale-relations are shown to be useful tools to obtain robust diagnostics, as examples: (a) HafH/3 ratio varies 
from the usually assumed value of 2.86, these variations imply the existence of a lower limit to the attainable precision in extinction 
estimations of AE(B — V) ~ 0.1.; (b) EW(Hfi) is a good age indicator with a small dependence on cluster mass, while £VK([Om] 
5007) shows a noteworthy mass dependence; (c) abundance estimations from ^23 are practically unaffected by variations of the cluster 
mass; (d) estimations from 523 and rf would improve if the cluster mass dependences were considered and (e) [On] 3727/Ha is a 
good star formation rate indicator for ages older than ~4.5 Ma. We also show that the ionizing cluster mass dependence explains why 
empirical calibrations produce more reliable diagnostics of some emission lines than photoionization models grids. Finally, we show 
preliminary results about the contribution of low mass clusters (M < 10 4 M Q ) to the integrated emission line spectra of galaxies, 
which can be as high as 80% for some relevant lines. 

Key words. HII regions - Galaxies: star clusters - Galaxies: stellar populations 



1. Introduction and motivation 

The intrinsic properties of star forming regions are of- 
ten inferred from observations of emission-line indicators. 
Examples of these indicators are the equivalent width of H/3 
(EW( Hfi)) for age estimation s, the radiation softness parame- 
ter 77/ dVflchez & P agel 1988) for cluster eff ective temperature 
(r e ff) esti mations, or the R23 (also ca lled O23) dPagel et al.lll979t) 
and S 23 dVTlchez & Esteban 1 19961) parameters as abundances 
indicators. So me reviews that deal with t hese indicators have 
been w ritten bv lFerlandl (120031) . IStasiriskal d2004 and lStasinskal 
(120071) . 

These and other indicators are us ed in estimation-met hods 
and/or classification methods , like the Baldw in et al.1 (11981b and 
Veilleux & Osterbrock ( 1987) diagnostic diagrams. These meth- 
ods make use of only a handful of emission lines, more specifi- 
cally, Ha, H/3 and some collisional lines of ions of O, N and S, 
so the understanding of the dependence of the intensities of such 
lines on the physical properties would allow us to calibrate the 
estimation methods. 

When using theoretical methods for such indicator- 
calibration purposes, photoionization models and stellar pop- 
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ulation synthesis models prove to be very useful tools, as 
they allow the control of relevant physical parameters and 
return the possible emission line spectrum associated with 
such physical conditions. Nevertheless, one of the main is- 
sues that limits their power is the difficulty of establishing 
whether the solution found is unique for a given observa- 
tion. One way to deal with such problem is building tailored 
models in which as many observational constraints as possi- 
ble are assumed: an example of this strategy is offered by 
ILuridiana. Peimbert. & Leithererf d 19991) : iLuridiana & Peimbertl 
(2001) and iLuridiana et al.l d2003bl) . where models that simul- 
taneously reproduce the intensity of several emission lines ob- 
served through different apertures are computed. A different 
approach is to compute grids or sequences of models. In this 
case fewer parameters are used but they span a larger range 
of values; this allow to apply the grids or sequences to di- 
verse objects. Some examples c an be found in the works o f 
| Huchraldl977h.lQlofssonl dl 989b. iGarcfa- Vargas & Dfazl d!994l). 
IStasihska & Izotovld2003b and lDopita et al.l d2006b), among oth- 



One of the parameters that can be considered in the model 
grids for star forming regions is the ionizing cluster mass. Yet, 
in spite of being a fundamental physical property, little work 
has been done regarding the ionizing cluster mass. Examples 
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can be found in Garcia- Vargas. Bressan. & Diazl (fl995allbh . 
Stasin ska & Leithereij (Il996h. IStasinska et alJ J200lh . and in a 
more indirect way in lDopita et al.l d2006bl) . nevertheless the ex- 
plicit dependence with the mass is not discussed in these works. 

In this work we center our attention explicitly on the influ- 
ence of the ionizing-cluster mass on the emission-line spectra. 
There are two motivations for this choice. The first one is the 
long term goal of investigating the accuracy of the modeling of 
the emission line spectra of ionizing clusters by means of evolu- 
tionary synthesis models, and the influence of sampling effects 
in the initial mass fu nction (IMF) on the emission line spectrum 
jCervino et al. 2003, Paper I). In the first paper of the series we 
investigated when sampling effects in the stellar populations of 
the cluster cannot be neglected in the modelling of the emission 
line spectrum, and we showed that cluster masses must be larger 
than 1 3 M to obtain realistic results for relevant emission lines, 
10 3 M being the cluster mass limit where a single star is as lu- 
minous as the w hole cluster (that is, the Lowe st Luminosity Limit 
(LLL), see alsolCeryino & Luridiana! 12004. for more details). 
Cervino & Luridiana ( 2004) propose to use the mean values ob- 
tained from synthesis models in clusters at least 10 times larger 
than the LLL to minimize IMF sampling effects. In this second 
paper we aim to obtain a reference scale (calibrated with the 
cluster mass) that provide a mean reference value that allows to 
evaluate quantitatively the dispersion due to the IMF sampling. 
As it is known, and we will see below, the intensity of hydro- 
gen recombination lines roughly scale with the number of ion- 
izing photons, and hence with the number of massive stars and 
the mass of the cluster. However, a similar scale relation has not 
been established for the case of forbidden lines. In the last paper 
of this series, we will use these scale relations to evaluate the 
variance due to IMF sampling in the emission line spectrum as a 
function of the cluster mass. 

The second motivation is to investigate the emission line 
spectrum of galaxies by the composition of different Hn regions. 
In general, the emission lines spectrum of galaxies has been 
modeled by means of single photoionization models where all 
the stars are concentrated in a single point (whatever the stel- 
lar density o f such hypothetic point) and with a single ionized 
gas nebula dStasinska et alJ 12001b ICharlot & Longhettil 1200 ll 
as examples). Only few a ttempts have been done to consider 
an ensemble of clusters dJamet & Morissell 120081) and to in- 
clude an initial cluster mass function (ICMF) in the simula- 
tions (D opita et al.ll2006ah . It can be claimed that the form of 
the ICMlQ is such that each logarithmic bin of cluster mass con- 
tributes a similar amount to the total flux and massive clusters 
will ther efore tend to wash o ut the stochastic effects of low-mass 
clusters (Dopita et al. 2006b). However, this statement implicitly 
assumes that the flux of all relevant emission lines, including 
collisional ones, scale linearly with the cluster mass. Checking 
the validity of such statement is the second objective of this pa- 
per. 

To study these questions we have computed a grid of pho- 
toionization models with ionizing clusters of different masses. 
In this paper we do not consider sampling effects but the scale- 
relations of the emission line intensities with the ionizing clus- 
ter total mass. We have included in our study some extreme 
(unrealistic) cases of cluster masses (e.g. clusters with masses 
lower than 10 3 M , affected by sampling effects, or clusters with 
masses larger than 10 7 M Q , which are not consistent with the 
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assumption of an instantaneous burst scenario). Such extremes 
allow us to explore the range where the scale relations are valid 
in a photoionization model, despite that they cannot represent 
real clusters, as well as to check if the emission line spectrum of 
a non-active galaxy can be modeled as a single nebula. 

In the next section the properties of this grid of photoion- 
ization models are described and the differences with grids from 
other authors are discussed. In section [3] the results of the mod- 
els are analysed and the luminosity functions are computed. In 
subsections 13.2.11 and 13.2.21 the reliability of these functions is 
checked comparing it with the results obtained directly from the 
models. The implications that arise from the analysis are dis- 
cussed in section HI and finally the conclusions of this work are 
summarized in section[5] 



2. Modeling strategy: defining the grid 

In the following we describe the input parameters assumed for 
the grid, and the differences with the grids elaborated by other 
authors. 



2. 1 . Continuum shape and intensity 

In this work we have used the spectral energ y distributions 
(SEDs) computed with the evolutionary code by Cervino et al. 
d2002al) . All the comp utations assume an instantaneous burst 
of star formation and a Sal peterl d!955l) IMF slope in the mass 
range 2 to 120 M . We have adopt ed the evolutio n ary tracks with 
standar d m ass-loss rates bylSchaller et alJ(11992l) : ISchaerer et al.l 
dl993alTbh ; ICharbonnel et al.l d 1993b. and th e following atmo- 
sphere models: ISchaerer & de Koterl d 19971) Cq Star for main- 
sequen ce hot stars mor e massive than 20 M , ISchmutz et al.l 
(1 19921) for WR stars, and lKuruczl d!991l) for the remaining stars. 
No X-ray emission due to the conversion of kinetic energy or the 
presence of young Supernova Remnants has been included. The 
models have been computed for the five metallicities available 
in the Geneva track dataset (Z=0.001, 0.004, 0.008, 0.020 and 
0.040) and cover the age range from 0.1 to 10 Ma with a time 
step of 0. 1 Ma. The resulting SEDs have been used as an input 
to the photoionization c ode Cloudy version 08.00, last described 
bv lFerlandet ail (119981) . 

The continuum intensity has been defined in terms of the 
amount of photons able to ionize hydrogen, Q(H°). We have 
used £2(H°) values that correspond to amounts of gas trans- 
formed into stars at the onset of the burst ranging from 1 to 
10 7 M . The mass step in the most relevant range ( 10 3 - 10 6 
M ) is 0.5 dex while for lower masses the mass step is 1 
dex. Obviously the 1 M case is just an academic case, since 
the lower mass limit of the IMF is 2 M . Moreover, clus- 
ters with masses below 10 3 M are also academic cases, since 
the (aver age) SEP used wou l d not be representative of real 
clusters dCervino et al.l 120031 ICervifio & Valls-Gabaudl 120031: 
ICervino & Lur idiana 2004, 2006). Whatever the case, the ex- 
tension down to fictitious 1 M clusters allows to explore the 
proportionality of the resulting emission line spectrum with the 
incident ionizing flux, independently of how "real" the incident 
flux is. 

This strategy differs from the one by iDopita etal . (2000), 
who use the SED of a zero-age stellar cluster to build a 
grid where the metallicity and the ionization parameter of the 
model nebulae are varied. It is important to note that using 
a zero-age SED is roughly equivalent to using a single stel- 
lar atmosphere model with T e ff between 45,000 and 50,000 K 
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(c.f.lMas-Hesse & Kunthll99Tl lG arcfa- Vargas. Bressan. & Dfazl 
Il995ah . with the exact value depending on metallicity. One of 
the conclusions of these authors, specifically that Hn regions 
can be modelled with zero-age SEDs, contrasts with the phi- 
losophy underlying evolutionary synthesis models, which are in 
their essence evolutionary. Additionally, were the SED of ob- 
served Hn regions described by synthesis models at zero-age, 
then all observed Hn regions should have similar H/3 equivalent 
widths (with values around 500 A), which is clearly in contra- 
diction with the observations. The fact that their grid reproduce 
the observed trends in some diagnostic diagrams is a nice ex- 
ample of the non-uniqueness of the possible physical conditions 
that produce a given emission line spectrum. 

A strategy more similar t o our s is the one by 
Garcia- Vargas, Bres san. & Dfazl dl995albl) . who present 
their results in terms of ages, metallicity and initial amount of 
gas transformed into stars. Our modelling differs from theirs 
in (a) the metal mixture, especially for what concerns the N/O 
ratio, which in our grid varies with the metallicity, and (b) the 
geometry, which is spherical in our case and pla ne-parallel 
geometry in Garcia- Varg as. Bressan. & Dfazl d!995al) . 

Our grid is si m ilar i n appro ach to the one by 

IStasinska & Leithererl d 19961) and IStasihska etail d2001l) . 
Differences are in this case the metal mixture, and the higher 
resolution in age and cluster size in our grid with respect to 
theirs. 

Finally, iDopita et alj d200 6b) used for the grid a parameteri- 
zation on metallicity, age and on a parameter that depends both 
on the ionizing cluster mass and on the pressure in the interstellar 
medium (ISM). In this sense the cluster mass is considered but 
its effects are mixed with the ones of the pressure. There are also 
differences in the metal mixture, in the age covering, which in 
our work is larger and has higher resolution, and in the covering 
factor , which in their work varies with the age. 

2.2. The chemical composition mesh 

In this work we assumed a du st- free nebula and the so- 
lar composition listed in Hazy dFerlandl 20081 table 14), 
which corresponds to the abundan ces by Grev esse & Anders 
d!989h with the extensions by iGrevesse & Noels! d 1993b 
except for C, N, O and He. The solar oxygen abun- 
dance has been set to O/H = 4.91 x 10~ 4 , following 
lAllende Prieto. Lambert. & Asplundld200ll). The solar C/O rati o 
has been set to 0.5 ( All ende Prieto. Lambert. & Asplu nd 2002), 
yielding a solar C abundance of 2.45 x 10~ 4 . For metallicities 
other than solar, the abundances have been linearly scaled with 
Z except for He and N. For the He abundance, we have used the 
following relation: 



by: 



In the case of nitrogen, we have assumed the N / O ratio given 



'-rr = 0.0786+ 18.4 



He 
~H 



(1) 



which has been obtained fro m the value of t he prim ordial 
helium abundance obtained by Luridiana et al. (2003b), Y p — 
0.2391, and the slope of the He v s. O relationsh ip quoted by 
Peimbert, Peimbert & Ruizl d2000l) . A7/AO = 3.5. lDopita eTatl 
(2006b) assume for the He abundance a steeper relation given 
by: 



log(g) = -1.5 farlog(§)<-4. 

= l.llog(g) + 2.9 forlog(g)>-4. 



The relati on has been obtained by a fit by eye of the data 
presented in iPilvugin. Thuan. & Vflchea d2003l). It shoul d be 
noted that, as shown by Pilv ugin. Thuan. & Vflchezl d2003l) . an 
exact N/O ratio does not exist, since its value changes from 
system to system even for a fixed metallicity. This dispersion 
could arise naturally if the chemical evolution of a galaxy 
were driven by small star forming regions, so that the sam- 
pling eff ects of the IMF in the c hemical evolution become 
relevant dWhite & Audouzd 119831: IWilmes & Koeppenl 119951: 
ICervino & Mollall2002t ICarigi & Hernandez! 120081) . Given this 
situation, we think it is more meaningful to fit the results by eye 
rather than using a formal method. 

We have also performed a consistency check on the as- 
sumed abundances of the stellar tracks. As explained in the co r- 
responding papers dSchaller et al.ll 19921 ISchaerer et al.lll993allbl 
ICharbonnel et al.lll993l) . the tracks assume a helium content by 
mass given by: 



AY 

Y - = Y„ + — Z, 
p AZ 



(3) 



— = 0.0737 + 52.5 '° 
H 



(£)■ 



with Y p = 0.24 and AY/AZ = 3 for al l metallicities except for the 
case Z=0.040 dSchaerer et al .111 993aL twice solar) where AY/AZ 
takes a value of 2.5. Additionally, the rela tive ratios of heavy 
eleme nts follow the ones obtained by the lAnders & Grevessa 
(198 91) solar mixture, so mewhat different from the values by 
IGrevesse & Andersl d 19891) quoted before, for all the metallic- 
ities. The absolute values have b een scaled wit h Z except for 
Z=0.004 dCharbon neTet al.lll993h and Z=0.040 dSchaerer et al.1 
fl993ah where IGrevesse! d!99 lFl has been used. Under this as- 
sumption, a N/O ratio of 0. 15 is assumed for all the metallicities 
except forZ=0.004 and Z=0.040, where it is 0.13. This is clearly 
inconsistent with the mixture we assume in the nebular gas, but, 
since we are not able to quantify the influence of a change of 
the N / O ratio on stellar evolution, and we prefer to adopt a neb- 
ular metal mixture more consistent with the observational data, 
there is no solution other than maintaining the inconsistency and 
hoping that it only produces minor effects. 

Finally, we have obtained the oxygen abundance from the 
different set of tracks. The resulting O/H values for each metal- 
licity, Z, are summarized in tab led] where we also sho w the gas 
O/H and Z values assumed in Stasiris ka et al.l d200ll) for com- 
parison. 

Since there is no feasible way of solving the inconsistency 
between the abundances of the evolutionary tracks and those 
adopted for the gas, present in our and similar works, we limit 
to call the attention to it and present our model grid in terms of 
Z, without expressing the metallicity as a function of the solar 
value. It should be kept in mind that a direct comparison with 
other works must be done with extreme caution. 



2.3. Density, geometry, filling factor and covering factors 

The assumptions on the electronic density, the geometry, the fill- 
ing factor and the covering factor define the properties of the 



(2) 



The reference paper cites lGrevesse et al ] dl99ll) instead of this one. 
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4.39 x 10- 5 




2.46 x 10- 5 


0.004 


1.77 x 10- 4 


1.70 x 10~ 4 


9.82 x 10 5 


0.008 


3.64 x 10~ 4 


3.40 x 10~ 4 


1.96 x 10" 4 


0.020 


9.75 x 10~ 4 


8.51 x 10~ 4 


4.91 x 10~ 4 


0.040 


2.12 x 10- 3 


1.70 x 10~ 3 


9.82 x 10- 4 



Tabl e 1. Value s of O/H implicit in the evolutionary tracks, in 
IStasinska et alJ d200ll) grid and in this work for a given Z. 



model nebulae. For the case of the density and geometry we 
have assumed rather conservative values: a hydrogen constant 
density of 100 cirT 3 distributed homogeneously over an expand- 
ing spherical geometry in an ionization bounded configuration. 
The assumed density of 100 cirT 3 is typical for Hn regions and 
it is lower than the critical densities of the species used in this 
work. The internal radius of the cloud, R m , has been fixed for all 
the models to 2 pc (logfl(cm) = 18.790426). 

The covering factor D./4n is the fraction of the solid an- 
gle around the ionization source covered by the gas nebula. 
There is a large am ount of observati onal evidence (see e.g . 
Beckman et al.ll2000t IZurita et al.|2002t ICastellanos et al.ll2002L 
Iglesias-Paramo & Munoz-Tu ndn 2002, among others) that a 
substantial fraction of the ionizing photons escape from the ion- 
ized nebulae. This finding may indicate either that the covering 
factor is smaller than one, or that the nebula is density bounded 
(or a mixture of both situations). Although the predictions of 
photoionization mode ls differ in the two situations (see, e.g. 
iLuridiana et al.ll2003bl for a full discussion), we have chosen to 
model the escape of ionizing photons adopting a covering factor 
and assuming the nebula ionization bounded, rather than assum- 
ing it density bounded. This choice has been motivated by sim- 
plicity considerations, since it offers the advantage of preserving 
the proportionality between the number of ionizing photons and 
the intensity of Balmer lines, in contrast to what happens in den- 
sity bounded nebulae. This means that we can test L(H/3) estima- 
tions by synthesis models based on Q(H°). In this work we have 
assumed a covering factor equal to 1, although estima ted values 
of the covering factor range between 0.2 and 0.8 (Zurit a et al.l 
120021) . A gain, our choice is based on simplicity considerations, 
since our objective is to obtain the shape of the luminosity func- 
tions of the emission lines intensities but not to calibrate it. 

For each age, metallicity and mass, the photoionization 
model has been computed until a minimum electron tempera- 
ture has been reached. This minimum electron temperature is 
4000K for Z=0.001, 0.004 and 0.008 and 200K for Z=0.020 and 
Z=0.040. Once a solution has been found, it has been used as 
input for another iteration, and so on, until a stable solution has 
been found. 

Additionally, we have obtained the continuum flux result- 
ing from Cloudy, obtained directly with the punch continuum 
command. However, this continuum includes the luminosity of 
the resulting emission lines. To exclude the emission line con- 
tribution, we have obtained it using the punch outward contin- 
uum comman d and then subtract it from the continuum flux (see 
Ferland 2008, for more details on the use of this command and 
the punch continuum one). The resulting continuum has been 
used later to obtain the equivalent widths of the corresponding 
emission lines. 

We want to note that the choices made here, especially those 
concerning the geometry, density, filling factor and covering fac- 
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Age (Ma) 

Fig. 1. Hy6 equivalent width versus age. Grey scale colors rep- 
resent metallicities. Thin solid line: 1 M . Thin dashed line: 10 
Mq. Thin dotted line: 100 M . Thin long-dashed line: 10 3 M G . 
Thick dash-dot line: 10 4 M . Thick dotted line: 10 5 M . Thick 
dashed line: 10 6 M . Thick solid line: 10 7 M . 



tors are merely convenient choices: we are not suggesting that all 
Hn regions can be modelled with this set of parameters. We are 
interested here only in obtaining a rough estimation of the rela- 
tion of intensity of relevant emission lines with the stellar cluster 
mass and how these dependencies vary from line to line. This 
goal is the driver in our choice of the grid which drastically dif- 
fers from how other authors define their grids, since they address 
a different problem from the one discussed here. 

3. Results of the grid models 

In this section we analize the results from the models. First we 
center our attention on the two hydrogen recombination lines Ha 
and HjS, specifically on their use for age and extinction estima- 
tions. Next we derive mass-luminosity scale-relations from the 
models and apply the same method to some relevant collisional 
lines. Finally we check the scale-relations results with the ones 
obtained directly from models. 

3.1. Hydrogen recombination lines 
3.1.1. H/5 equivalent width 

The resulting EW(H/3) from the grid is shown in Fig. Q] The 
plot shows the results for the complete grid (including all the 
assumed cluster masses). Different line types account for clus- 
ter masses while grey scale colors represent metallicities. For a 
given metallicity, the curves at different masses are similar with 
just slight differences that depend on the ionizing cluster mass. 
Excluding the 1 M cases, these differences are always less than 
0. 1 dex. The 1 M cases diverge from the rest of models at old 
ages because the computation goes deeper into neutral gas and 
the continuum at H/3 wavelenght increases. 
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Fig. 2. Ratio of the HB luminosity to the number of hydrogen 
ionizing photons, Q(H°), versus age for models with M > 10 3 
M . The shaded zone corresponds to the range of values of the 
ratio. Dash-dot-dot line: value assumed by synthesis models. 



As a safety check to the resulting EW(H/3) from synthesis 
models computations, we have compared the ratio of the HB lu- 
minosity with the number of hydrogen ionizing photons, Q(H°). 
Synthesis models assume a putative constant value around 
L(Hff) = 4.8 X 10- 13 (Q/47r)Q(H° ) (e.g. ICervino & Mas-Hessd 
ll994tlLeitherer & Heckman Il995h for T e = 10 4 K, n e = 100 cm" 3 
and He/H = 0.1. Figure [2] shows that the resulting conversion 
value for M > 10 3 M only varies ~20% during the considered 
ages. 

Summing up, EW( HB) is a "good eno ugh" indicator of the 
age of an Hn region dCopetti et alj 119861) once the covering 
factor is known. EW(HB) d epends mainly on the metallicity 
( ICervino & Mas-Hessd fl994) and age of the cluster with almost 
no variations due to the size of the cluster, if sampl ing effects are 
not considered (but see ICervino et al . 2()()(). l2()()2hl) . 



3.1 .2. Ha, HB and extinction 

More interesting is the behaviour of the HafHB ratio, used com- 
monly as an extinction index, and for which a constant ratio of 
2.86 is normally assumed. Figure [3] shows such ratio, obtained 
for models with M > 1O 3 M , the line types have the same mean- 
ing as in Fig.Q] The ratio computed for case B and T e = 10 4 K is 
also shown for comparison. The figure demonstrates that in fact 
such an assumption is incorrect. The largest discrepancies be- 
tween the results of the models and the value normally assumed 
occur generally for Z=0.040. Nevertheless, the largest ones cor- 
respond to models with Z=0.008, 1 M Q (not shown in the figure) 
and ages around 10 Ma, which are ~50% greater than the usual 
assumed value and lead to E(B — V) = 0.36. For that cluster 
mass the mean spectrum obtained by synthesis models is just 
an academic result, as we pointed out before. Considering only 
the models with M > 1O 3 M a maximum discrepancy of 25.8% 
with respect to the ratio normally assumed is found for models 



Fig. 3. Ha/HB ratio versus age. Dash-dot-dot line: normally as- 
sumed value. Rest of symbols as in Fig.Q] 

with Z=0.04, 10 Ma, and masses larger than 1O 6 M . However, 
this only sets an upper limit to the bias since at such ages the 
Hn region are so faint that they are not detected. Considering 
lower and more relevant ages (0.1-6 Ma), differences near 18% 
are found for Z=0.04 and masses 10 3 5 and 10 4 M at 6 Ma. If 
we also exclude the highest metallicity models, differences near 
9% are found for Z=0.02, M > 4.5M at ages near 6 Ma. 

The use of an incorrect Ha/HB ratio leads to an incorrect 
extinction estimation which affects the line intensities inferred. 
For example, a difference of 9% with respect to the normally 
assumed ratio leads to a spurious E(B - V) - 0.08. Although 
this value is small, this exercise shows that E(B - V) estimations 
from the emission line analysis (without the use of tailored mod- 
els) suffer from intrinsic uncertainties which may be as large as 
0.1. This effect is independent of the extinction due to the Hn 
region dust content. If dust effects were taken in account, then 
the scatter in the HafHB ratio would be different: bigger Hn re- 
gions would have a larger amount of dust, the dependence with 
the geometry would be stronger and a tailored analysis would be 
needed. 

3.1.3. Scale-relations 

In order to study the variation of the intensities of the emission 
lines considered in this work, we made ordinary least square lin- 
ear regressions of the logarithm of the line intensities versus the 
logarithm of the cluster mass for each line, metallicity and age. 
The linear fits provide us scale-relations between the line inten- 
sities and the cluster mass of the form 



log(L) = a +y Slog(M), 
or equivalently 

L = AxM /3 , 



(4) 



(5) 



where A=10 ff . These scale-relations, as said above, have been 
built to obtain a first estimation of the behaviour of the emission 
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Fig. 4. From top to bottom: slopes, log(L[erg s M ']) and 
correlation coefficient of linear regression of log(L(H/?)) versus 
log(M). Grey scale colors indicate metallicities as in Fig.Q] Dots 
indicate models with age equal to 0.1, 1, 2, 3, . . . Ma. 

lines with the mass of the cluster. Due to the assumptions made 
in the models and the simplicity of the description with power- 
laws, we strongly discourage from using this method when fine 
results are sought. Let us stress again that our goal is just to ob- 
tain an approximate quantitative representation of the H n region 
luminosity function. 

The regressions have been computed only with models with 
M > 10 3 M o . This condition to th e cluster-mass agrees with the 
minimum cluster mass obtained bv lCervino et al . (2003) for a re- 
liable (although affected by sampling effects) description of the 
ionizing continuum by the mean ionizing continuum obtained in 
synthesis models. Regressions with models in all the mass range, 
including the purely academic case of clusters less massive than 
1O 3 M , have also been computed but they have been discarded, 
because the agreement in diagnostic diagrams between regres- 
sions and models in the high mass range is worse than when 
only models with M > 1O 3 M are used (see section [J.2.U . 

We have also studied an a posteriori additional constraint 
on the observability of emission lines. This constraint has been 
parameterized in terms of a distinction between lines with EW 
smaller and larger than 1 A respectively. Of course, this is an 
optimistic value: if older stellar populations were present, the 
continuum level would be higher and the resulting EWs lower. 
However, we assume just a single stellar population in an in- 
stantaneous burst mode of star formation in our modeling. The 
effects of these additional populations on the continuum level 
and EWs measurements can be modeled through the use of a 
different star formation history, which goes beyond the scope of 
this paper. 

The EW computations have been carried out using the 
line luminosity and the continuum luminosity as computed by 



Cloudy, which takes into account the two contributions to the 
continuum luminosity (stellar and nebular continua). We have 
obtained continuum scale relations using Eq. [4] As expected, 
the continuum luminosity (including the nebular contribution) 
scales linearly with the cluster mass, although there is a small 
dispersion when the continuum obtained from the scale relations 
is compared with the continuum from particular models (around 
7%, compatible with the scatter in the hydrogen recombination 
lines shown previously). We have computed the EW of the emis- 
sion line using the scale relations following: 

EW ^ AjteXMP = Aline xAf8 -l | (6) 

A C ont x M A cont 

where Aij ne , following Eq. [3] is the luminosity of the line when 
the cluster mass is 1 M Q , obtained from the extrapolation of the 
power law down to this value. Analogously, A cont is the luminos- 
ity of the continuum when the cluster mass is 1 M Q . Since we 
are interested in the relation of the luminosities with the cluster 
mass, we center our attention on the slopes of the scale-relations. 
However, the EW obtained from the scale-relations are more de- 
pendent on the geometrical assumptions (inner radius, covering 
and filling factors) than the scale-relations slopes themselves, 
due to the large dependences of Ai; ne and A cont on the geometry. 
So, we will use this constraint on EW as a self-consistency check 
by comparing the results from models that fulfill EW > 1 with 
scale-relations results with EW estimations obtained from the 
scale-relations themselves. The values of the normalized lumi- 
nosities will be shown in a posterior work when a proper ICMF 
and star formation history will be included for models of non- 
active emission line galaxies. 

In Fig. [4] the slopes, the luminosity normalized to the total 
mass and the correlation factor (plotted as 1 - r) as a function 
of the age obtained for H/3 line are shown. The slopes for H/3 
intensity are practically constant and with a value almost equal 
to one. This behaviour is expected for hydrogen recombination 
lines (the Ha intensity shows similar results) because their in- 
tensity is practically proportional to the 2(H°) of the ionizing 
cluster, which is directly related with the cluster mass. This is 
coherent with the small range of the L[H/3)/Q(H°) ratio shown in 
Fig.|2]for each age. Also, 1 — r « accounts for this proportion- 
ality. 

Lines that have power-laws with /3 ~ 1 are particularly in- 
teresting because, as we have said, they vary linearly with the 
cluster mass, or equivalently with Q(H°), just as the continuum 
produced by the stellar cluster. In general, the use of line ra- 
tios with lines that have similar f3 values provide indices that are 
independent on the cluster mass, and so, such indices can be cal- 
ibrated with observations of clusters without the need to worry 
about the cluster size (if IMF sampling effects are not taken into 
account). This is the case of EW(H/3), as can be seen in Fig. 
Q] because the intensity of the continuum is proportional to the 
cluster mass. 

3.2. Collisional lines 

Following the same method described before we have obtained 
scale-relations for some relevant collisional lines. In Fig.[5]the 
slopes for [On] A 3727 A, [Ora] A 5007 A, [Nn] A 6584 A, 
[Sn] A 6716 A, and [S m] A 9069 A are shown. The grey scale 
colors of the lines have the same meaning as in Fig.Q] We have 
represented with dotted lines the ages where the scale relations 
computed through Eq.|6]do not produce EW > 1 A for all masses 
in the 10 3 - 10 6 M mass range (we have excluded the 10 7 M Q 
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Fig. 7. [O in] 5007/H/? vs [N n] 6584/Hff diagnostic diagrams for each metallicity. Grey scale lines represent results from models and 
scale-relations for different cluster mass. Narrow lines with symbols represent the age evolution of models that fulfill the EW > 1 
A criteria; different symbols correspond to different ages. Dotted/solid lines are the results of applying the scale relations to the 
individual cluster masses, being the dotted part the region where the EW (as computed by Eq. |6]l of any of the lines involved is 
lower than 1 A. For reference we also show Ivan Zee et al.l d 1998b observations (open black squares), and the iKewlev et all (1200 ll) 
curve (long dashed) with its error limits (short dashed curves). 



case, since it is a purely academic case with no implications for 
observations). For these emission lines r > 0.98 in all the age 
range, therefore within our models the logarithm of the intensity 
of the emission lines is, to a good approximation, linear with the 
logarithm of the cluster mass. 

The slopes of [On] A 3727 A, [Nn] A 6584 A, and [Sm] A 
9069 A have a similar behaviour. In the three cases, there are 



two regimes, a young one with slopes around 0.7-0.8, and an a 
older one with slopes around 1 (excluding the Z=0.04 models). 
The transition between the two regimes happens in the age range 
4 - 5.5 Ma and it is due to a softening of the ionizing continua 
in such age rang e. Comparing with the evolu tion of the cluster 
T e fi, defined as in Mas-Hesse & Kunfhl(ll991h . the change in the 
slope is associated with a cluster T e g ~ 3.5 x 10 4 K following 
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Fig. 8. Om 5007/H/3 Vs S n 67 16,673 1/Ha diagnostic diagrams for each metallicity. Symbols as in Fig. [7] 



the results from|C ervino & Mas-Hessel d 19941) and lCervino et alj 
( 2002b Such cluster T e ff is not due to the predominance of a 
specific type of star but to a mix of different types of stars. 

Since the different values of the slopes indicate different de- 
pendences with the mass, and taking into account that the contin- 
uum is proportional to the cluster mass, it is straightforward that 
the equivalent widths also present different dependences with the 
mass. This can be seen in Fig. [6^ in which £W([On] 3727) is 
shown for M > 1O 3 M and all the metallicities. There is a clear 
dependence of EW{ [On] A 3727) with the mass for young ages 



3 http : //www . laef f . inta . es/users/mcs/SED/ index . html 



(slopes different from 1). On the other hand, the EW is almost 
independent on the cluster mass for old ages. 

[On] A 3727 A and EW( 3 727) are used a s star 
formation rate (SFR) indicators (jGallagher et all j 1989; 
Kennicutt [1991: I osa-Gonzalez et all u002, among oth- 
ers). The method is based on linear relations between 
[On] A 3727 A and Ha or H/3. Many authors have stud- 
ied the dependence of these relations on the luminosity 

■ h I I I ■ I I J 

ansen et al. 2001; Araaon-Salamanca et al. 2003), the red- 
denin g dRosa-Gonzalez et al.l |2002[ Aragon-Salamanca et al j 
20031) and abundance s dJansen et ^1.1 1200 lHKewlev et al l 12004; 
Mouhcine et all 120051: iMoustakas et al.ll2006l) . From the value 
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of the slopes for [O n] A 3727 A and Ha (being the last similar 
to the ones shown in Fig. 0), we obtained that the [O n] 3727/Ha- 
ratio roughly depends on M~ 3 in the 0.1-4.5 Ma range, having 
a dependence flatter than M° 05 for larger ages. Therefore at very 
early ages the SFR estimation has a non negligible dependence 
on the mass of the ionizing cluster of the observed Hn regions. 
However, it is a good SFR indicator for older ages. Our study 
suggests that a recalibration of [On] A 3727 A using only 
old (age > 4.5 Ma) clusters would improve its use as an SFR 
indicator and that the method should be applied to galaxies 
dominated by H n regions older than 4.5 Ma. 

It is noteworthy that the slope of [O in] A 5007 A is nearly 
1. Considering all the metallicities except Z=0.04 and ages 
lower than 7 Ma, we can adopt a value of 1.05 as a repre- 
sentative slope for this emission line. In such case we obtain 
EW([Om] 5007) oc M 005 . Th erefore EW(\Om\ 5007) seem s to 
be a suitable age indicator as IStasihska & Leithererl d 19961) had 
proposed. Nevertheless the dependence with the cluster mass is 
large enough to produce a non-negligible uncertainty in the age. 
This can be checked in Fig. [6J3 where the equivalent width of 
[Om] A 5007 A is shown. Comparing it with Fig. [1] it can 
be seen that the sequences for equal metallicity span a larger 
£W-range than for EW(H{3). Moreover there are more situations 
of age degeneracy if [O in] A 5007 A is used. Therefore cau- 
tion is needed when using this line for age estimations. In Fig. 
|6j? it can also be seen that the age range for the application of 
EW([0 in] 5007) as an age indicator is smaller than for EW(H/3), 
because EW([0 m] 5007), and therefore the detection of the line, 
decreases faster than EW(Hfi). 

The slopes of [S n] A 6716 A are smaller than 1 except for 
Z=0.04 at older ages. This means that the EW and ratios of [S n] 
A 6716 A with Ha and Hp will decrease with increasing cluster 
mass. 

The reliability of these scale-relations can be checked com- 
paring the results obtained with them with the results obtained 
directly from the models. The comparisons have been made with 
some diagnostic diagrams and computing some relevant param- 
eters. At the same time, clues of the influence of the cluster mass 
on such diagrams and parameters have been obtained. 

3.2.1. Diagnostic diagrams 

We have used two diagnostic diagrams to compare the scale- 
relations with the models. Let us remind that scale-relations built 
with models in the whole mass range (1-10 7 M ) do not agree 
well with the results obtained directly from the models, espe- 
cially in the high-mass range where the luminosity is overesti- 
mated. The inclusion of low mass clusters, which is just an aca- 
demic case, produces a concave curvature in the logM - logL 
plane, so the inclusion of low mass clusters produces a steeper 
slope (and a lower normalized luminosity). As we explain in 
section I3.1.3I we primarily seek a good agreement in the high- 
mass range. For this reason our final scale-relations have been 
computed only with models with M > 10 3 M . This clearly im- 
proves the results of scale-relations in the high-mass range; fur- 
thermore, we have also checked that for low masses the results 
are acceptable (though overestimated). In the following, we fo- 
cus only on the simulations with clusters larger than 10 3 M . 

In Fi g. [Tithe [O ml 5007/H /3 vs [N n] 6584/Ha diagram pro- 
posed by Bald win et al.ldl981l) is shown. Narrow lines with sym- 
bols represent the age evolution of models that fulfill the E W > 1 
A criterion. Different symbols correspond to different ages. Note 
that the models do not fulfill the EW for all ages (and they have 



not been plotted in such cases). Dotted/solid lines are the re- 
sults of applying of the estimations of the line intensities with 
the scale relations obtained previously to the individual cluster 
masses, the dotted part being the ages where the EW of any of 
the lines involved in the diagram is lower than 1 A. The t heoret- 
ical upper limit for starbursts from iKewlev et alJ (l200lh (black 
dashed line) and its estimated errors (black dotted line) are also 
shown for reference. 

In general, there is a good agreement between the models 
and the results obtained using the scale relations, and there are 
only noticeable discrepancies for the Z = 0.04 case. It is also in- 
teresting to find in the Z = 0.001 case that only the less massive 
clusters produce £W([Nn]) larger than 1. It is consistent with 
the results of Fig. [5] which show that, given the slope lower than 
1 of the L([Nn]) scale relation, the /sW^Nn]) decreases as the 
cluster mass increases for young ages. Our models and scale re- 
lations show that for fixed metallicity and T e g, i.e. constant age, 
there is a sequence in decreasing Q(H°), or equivalently in de- 
creasing mass. In general, models and scale relations have trends 
similar to, or at least compatible with, the observations. 

Some of our s c ale-r elations (and models) are slightly above 
the IKewlev et al.1 (l200ll) curve. Such models correspond to 
clusters of 10 7 M and metallicities in the 0.004-0.02 range. 
IKewlev et al.l (12001 1) estimated an error of +0. 1 dex due to po- 
tential errors in the modelling. Taking into accou nt such errors, 
all ou r scale-relation estimations agree with the IKewlev et al.l 
(120011) limits. 

In Fig.|8]the [Om] 5007/H/3 vs [S11] 6716+6731/Ha diag- 
nostic diagram is shown. The results estimated from the power- 
law scale-relations follow the trends of the models, and also re- 
produce some of the features (peaks, turn-off and EW) that the 
sequences of models present. The diagram shows clearly the ef- 
fects of the cluster mass in the [S n] 6716+6731, which produce 
a serie of parallel sequences. Again, there is a good agreement 
between models and scale relations, with the larger discrepan- 
cies at Z = 0.04. Scale relations show that it is more difficult to 
detect the [S 11] lines in more massive clusters. Since the slope 
in the log M - log L relation is almost constant for all ages and 
lower than one, both the [S11] 6716+6731/Ha and the EWofthe 
line decreases with increasing cluster mass. In this diagram, ob- 
servations and models/scale relations have a poorer agreement 
than in the [O m] 5007/H/3 vs [N 11] 6584/Ha- diagram. However, 
notice that [S 11] AA 67 1 6, 673 1 A, as [0 1] A 6300 A, are likely to 
be affected by shocks and since Cloudy does not include shocks, 
the results for those lines are not complete. 

From these comparisions we conclude that, although not 
perfect, the scale-relations are suitable tools for cluster masses 
larger than 10 3 M . 

3.2.2. The 77 parameter, R 2 i and S23 

The 77 parameter was defined by Vflch ez & Pagell dl988l) as: 

o + /o ++ 



' S+/S++ ' 

and it is related with the observational ratio: 

, _ [On] 3727 + 3729/ [Om] 4959 + 5007 
77 ~ [S 11] 6717 + 6731/[S in] 9069 + 9532 ' 



(7) 



(8) 



The slopes of the scale-relations computed for each of the 
lines involved in the ratio can be used to estimate the 77' depen- 
dence on mass. As can be seen in Fig.|5j for each line the com- 
puted slopes have similar values and trends for all metallicities 
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,1.05. 



except for Z=0. 04. Thus, as a rough approximation, we can de- I{[0 iii]5007) oc M 
fine for each line a representative slope and obtain the following 7([Sn]6716) oc M ' 73- 
proportionalities valid for all the metallicities considered except ( 77 



forZ=0.04: 



I([S m]9069) oc 



M 



1.03 



if age < 4.5 Ma, 
if age > 4.5 Ma. 



(10) 
(ID 

(12) 



. M 015 if age<4.5Ma, 
7([Oll]3727) oc < Ml03 . f age>45Ma . 



(9) 
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[O in] A 4959 A, [S n] A 6731 A and [S m] A 9532 A present 
dependences with the cluster mass similar to their associated 
lines from above. Using these proportionalities we obtain for 77': 



AT 
Jvt 



0.2S 



if age < 4.5 Ma. 
if age > 4.5 Ma. 



(13) 



These trends can be seen in Fig. [9] where 77' versus log(Z) 
is represented for each cluster mass. The grey scale colors, un- 
like in previous figures, indicate the age. Results obtained di- 
rectly from models are shown with circles connected by narrow 
lines. We have only plotted the models that fulfil the EW criteria. 
Results computed with scale-relations are represented with dot- 
ted/solid lines with a larger line weight. Dotted parts of the lines 
connect metallicity points where either point does not fulfill the 
EW criteria as obtained by the scale laws. The agrement within 
the models and scale relations is remarkable. At ages lower than 
4.5 Ma, scale-relations produce 77/ values that decrease with in- 
creasing mass, while at larger ages 77/ increases with mass. This 
behaviour is due to the variation of the scale relation of the [S m] 
and [On] lines. 

The R23 paramete r was proposed as an abundance calibrator 
by Pagel etal] d!979l) . and is defined as: 



^23 



[0 11] 3727 + [O hi] 4959 + 5007 
H/3 • 



(14) 



Following the same method as for 77', we estimate the de- 
pendence of R23 with the cluster mass. Taking into account that 
H/3 oc M and using Eqns.l9land[T0l we obtain: 



R 2 3 <* 



l/?23 



M 



-0.25 



0.05 



A R23 M U 



-B*23M° 05 



if age < 4.5 Ma. 
if age > 4.5 Ma. 



(15) 



where the A«23 and B«23 factors are obtained from the intensity 
of the lines used in R23 normalized to the total mass, therefore 
they depend on age and metallicity. 

Recalling that the continuum has a slope equal to 1 (see sec. 
13 - 1 - 3b the same result is obtained if the EW of the lines are used 



Kobulnicky & Phillips (2003) proposed. 



In Fig. [10] the values obtained for ^23 are shown. The de- 
pendence in mass estimated above, in spite of being small, can 
be appreciated when comparing carefully the plots for different 
masses. The agreement between scale-relations and models is 
better than for the 77' parameter. 

The S 23 parameter was defined bv lVflchez & Estebanl (fl996) 
in a way similar to R23: 



>23 



[S n]6717 + 6731 + [S m]9069 + 9532 



(16) 



For this parameter we obtain a mass dependence of the form: 



>23 



A s2 3M-°- 21 + B S 23M-°- 23 



A s2 3M-°- 21 + B S 23M' 



0.03 



if age < 4.5 Ma. 
if age > 4.5 Ma. 



(17) 



where the As 23 and B.S23 factors depend on age and metallicity. 

As R23, S 23 also shows a dispersion of the results of the mod- 
els due to different cluster masses (see Fig.fTTTi. A relevant fact 
is the variation of S 23 with the age through the mass range. As 
the mass increase the sequences for older ages rise and the lower 
ages ones go down. This behaviour can be explained with the 
different dependences of S 23 with the mass for ages below and 
over 4.5 Ma. 



4. Summary of results and discussion 

Using models of Hn region ionized by clusters of different 
masses, we have studied the variation of the H 11 region emis- 
sion with the ionizing cluster mass. For this task, scale-relations 
in form of power-laws between the intensities of some of the 
most relevant emission lines and the ionizing cluster mass have 
been obtained. The reliability of these power-law scale-relations 
has been checked by means of diagnostic diagrams. We check 
that the general agreement between power-law scale-relations 
and models is better if we compute the power-laws using only 
the models with M > 10 3 M , with the intensities of collisional 
lines being overestimated by the power law relations in the case 
of clusters with masses lower than 10 3 M . However, the ion- 
izing flux of such low mass clusters cannot be described by the 
average ionizing flux obtained by synthesis models (c.f. Paper 
I). 

We want to stress that, due to the assumptions in the mod- 
elling and the simple linear approximation, these scale-relations 
are just a partial first approach from which we can obtain an 
overall view of the dependence of line intensities with the cluster 
mass, but these scale relations must not be used as an exact rep- 
resentation of the variation of the intensities of the emission lines 
with the cluster mass. In this sense, our scale-relations, although 
incomplete, are very useful tools that describe mean luminosity 
functions. Therefore they could be used when considering the 
mass dependence of the properties of massive H 11 regions sam- 
ples, but should not be applied to individual objects. An example 
of the limitations of our scale-relations is that they cannot predict 
the variation of the Ha/HB variation of the models with respect 
to the normally assumed value, variation that induce a spurious 
extinction effect. 

In spite of these limitations, simply estimating the trend of 
the luminosity distribution function allows to make interesting 
predictions. The exponents of the power-law scale-relations in- 
dicate us the mass dependence of the line luminosity. With this 
dependence it is possible to estimate the influence of the clus- 
ter mass on EWs and on some relevant line intensities ratios. 
If a quantity has no dependence on cluster mass it can be cali- 
brated with observations with no need of age and cluster mass 
estimations, otherwise a correction in mass should be applied. 
We check that EW(HB) is almost insensitive to mass, while 
^W^Oiii] 5007) has a small dependence. This implies that 
ZsWIOiii] 5007) is not as good as an age indicator as EW(HB), 
and that we have to take into account the cluster mass when esti- 
mating ages with EW{[0 111] 5007). We also explore the variation 
of fiW^On] 3727) with mass cluster and we obtain an impor- 
tant mass dependence for ages lower than ~4.5 Ma, depending 
on metallicity, and almost no mass dependence at older ages. 
The different mass dependences of [O n] A 3727 A, [O in] A 5007 
A, and therefore of their EWs, offer a possible explanation for 
the different disper sions in their respective log(£ 'W r /£'W(H/3)) vs 
log(I/T(H/3)) plots (Kobulnic kv~& Philhpsll2003l) . The mass de- 
pendence of the [On] 3727/Hor ratio, similar to £W([On] 3727), 
implies that the [0 11] A 3727 A intensity is not a good SFR in- 
dicator if clusters with ages lower than 4.5 Ma are present. This 
suggests that the [0 11] A 3727 A intensity should be recalibrated 
and that its use should be limited to galaxies dominated by H 11 
regions older than 4.5 Ma. 

The power-law scale-relations also allow us to identify emis- 
sion line ratios with small mass dependence. This is the case of 
[On] 3727/[Nn] 6584 rati o, pr oposed as an abund ance indica- 
tor bv lvan Zee etal] d!998l) and lDopita etal] d2000t) . The values 
and behaviour of the slopes for these two lines are similar, with 
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greater differences at the highest metallicities. Therefore their 
ratio will have a small mass dependence at any age. This can be 
seen in Fig. [12] where the [Nn] 6584/[On] 3727 ratio is plotted 
versus the age. The thin sequences obtained for low metallici- 
ties account for the small mass dependence. This clearly helps 
to separate the sequences for differents metallicities, allowing a 
better abundance estimation. The drop of [Nn] 6584/[On] 3727 
for Z=0.040 (long-dashed red line) at ages older than 6 Ma is 



due to a higher contribution of recombination emission to the 
[On] A 3727 A luminosity. 

The scale-relations have been used together with the defini- 
tion of some parameters (77/, ^23 and S 23) in order to estimate 
their dependence with the cluster mass. For 77' and S 23 the mod- 
els show variations of these ratios with the cluster mass that, in 
spite of the strong assumptions in the modelling, stress the need 
of taking into account the mass of the ionizing clusters when 
estimating physical properties from these ratios. For M > 10 3 
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Fig. 11. S23 for each cluster mass. Symbols as in Fig. [9] 



M , 77' and S 23 have opposite mass dependences at ages earlier 
and older than 4.5 Ma. On the other hand, R23 shows almost 
no dependence on mass. In this sense, R23 seems a quite robust 
indicator. Moreover, it also means that R23 could be calibrated 
without having to worry about the mass of the ionizing clusters, 
and therefore objects of different masses can be used together in 
its calibration. 

Additionally, the scale relations allow us to estimate when an 
emission line can be detected taking into account the continuum 



level. It turns out that it is easier to detect emission lines in low 
mass clusters when the scale relations have slopes lower than 
one. 

It can be argued that the variation with the cluster mass can 
also be described by variations of other parameters, but the use of 
the mass as a parameter allow to use directly the synthesis mod- 
els and an easier approach to other properties parameterized with 
the mass (as the star formation history) and to the the IMF sam- 
pling problem. Moreover, implementing the mass dependence 
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in the calibration of an estimation method would improve the 
estimations. Since including the cluster mass in H n region mod- 
els is a non-trivial task, and it is not a sufficiently studied issue, 
the theoretical calibrations of estimation methods with mass de- 
pendent indicators should be used with extra caution. Therefore 
empirical or semi-empirical calibrations would work better than 
theoretical ones since the former consider in a natural way the 
variation in the cluster mass, although it is necessary to estimate 
the mass by alternative methods (e.g. star counts or dynamic ef- 
fects). 

As we have said, our models and scale-relations are less re- 
alistic for the lowest masses (M < 10 3 M ). To improve the 
reliability in general, and in the low mass range in particular, 
an estimation of th e influence of sampling effects of the stel- 
lar IMF is needed (ICervino & Mas-Hesselll994t ICervino et alJ 
2000, 12002bl 120031 ICervino & Luridiandl2006Tr This would al- 
low to describe properly the possible distribution of the cluster 
ionizing continua. Such distribution for low mass clusters has a 
large dispersion and bimodality in T s s d Villaverde et al.l 12009). 
That imply that there could be low mass clusters with an ioniz- 
ing flux equal or greater than a larger mass cluster. The effects 
of the dispersion of T e g for low mass cluste rs on the H n region 
spectra are treated on lVillaverde et al](l2009l) to which the reader 
is referred. 

We can obtain the same back-of-the-envelope estimation of 
the influence of low mass clusters on the emission line spec- 
tra of a non- active galaxy, and obtain a first order check of the 
statement by iDopita et al.l (l2006bl) on the lack of influence of 
low mass clusters (and h ence, sampling effec ts) on the total line 
flux of collisional lines. Domta etal.1 <|2006a) assume a lognor- 
mal ICMF 



P(M) 



1 



Ma- V27T 



2<r2 



(18) 



with <x = 1.7 and /j — In 100 M Q . Such a distribution has a mean 
about 424 M Q , a median equal to 100 M and a mode (maxi- 
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Fig. 13. Contribution of clusters less massive than 10 4 M to the 
total luminosity as a function of the slope of the mass luminosity 
relation for different ICMF recipes. 



mum value) about 5 M . However, in Dopita et al. (2006a) they 
argued than they only make use of this distribution for cluster 
masses larger than 10 3 M and that massive cluster wash out the 
stochastic effects. 

We test their statement by computing the fraction of the to- 
tal luminosity of the emission line that is produced by clusters 
with mass less than M x = 10 4 M . The use of 10 4 M instead 
of 10 3 M is because for 10 3 M clusters, synthesis models in- 
cludes stars more lu minous than the cluster itself (the Lowest 
Lumi nosity Limit, see ICervino et al J2003t lC ervifio & L uridiana 
2004, for details), so at such masses there is not only a 
sampling effect problem, but a physical one in the modeling. 
Clusters in the 10 3 - 10 4 M are still affected by sampling ef- 
fects with luminosity d istributions that are not gaussian (c.f. 
ICervino & Luridiana 200^), so we take 10 4 M as a conservative 
lower limit for clusters affected by sampling. The integration of 
a lognormal ICMF with a power-law mass dependent luminosity 
is: 



AM a P(M)dM = --e ofI+ °- 



Erf 



\i + acr 



InM 



V2cr 



(19) 



with Erf(x) the gaussian error function. Hence, the contribution 
of clusters less massive than M x as a function of the slope of the 
mass-luminosity relation is: 
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We have also computed the same quantity for the case of a 
power-law ICMF, 



P(M) = CciM""", 



(21) 
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which produces the following contribution of clusters less mas- 
sive than M x as a function of the slope of the mass-luminosity 
relation: 

L( M < M x ) <T ad - M x +a - a « 



Figure Q~3] shows the variation of the contribution of clusters 
less massive than M x to the total luminosity as a function of 
the slope a of the mass-luminosity relation. The approximative 
values of a for different lines has been drawn as vertical lines. 
We have also assumed a M ]ow = 10 3 M and a M up = 10 6 , 10 7 
M as limit s of the ICMF. In the case of a lognormal ICMF as 
proposed by Dopit a~et~aT1 d200 6b). the figure shows that clusters 
in the 10 3 — 10 4 M range are responsible for more than 70% of 
the luminosity of the [O m] emission lines for all ages (and [O n], 
[S in] and [N n] at ages larger than 4.5 Ma). This contribution is 
larger than 80% for [S n] in all the age range (and [On], [S m] 
and [Nn] at ages lower than 4.5 Ma). This result is independent 
of the choice of the M up : the high mass tail of the lognormal 
distribution is steeper than the one of power law distributions, 
so high mas s clusters have a low impact on the total flux, in 
opposition to Dopit aet al.l d2006al) as sumption. In this situation, 
the statement by Dop ita et all d2006al) is clearly wrong. 

We have also show the case of a ICMF described by a power 
law with exponent a c \ = 2. In this situation, the contribution 
of massive clusters to the total luminosity is more important 
with a significant dependence of the choice of M up in the ICMF. 
However, the contribution of low mass clusters to some of the 
emission lines ([On], [Sin] and [Nn] at ages smaller than 4.5 
Ma) is around 50%, and only around 20% for [O in] for all ages 
and [O n], [S m] and [N n] at ages larger than 4.5 Ma. 

Although we have not taken into account either the star for- 
mation history or clusters less massive than 10 3 M in this calcu- 
lations, it is enough to prove that the statement about t he washing 
out of t he stochastic effects of low mass cluster by Dopit a et~al] 
(2006b) needs to be revised in detail. This subject will be ad- 
dressed in another work. 

Finally, we have compared our scale-relations with the 
Stasi nska et al.ld2001l) models. Those models considered several 
different values for density and filling factors, moreover the inner 
radius varies from model to model. For a consistence check we 
have only considered instantaneous burst models with IMF up- 
per stellar mass similar to ours. In Figs[l4][l5]and[T6]the intensi- 
ties of some of the most relevant emission lines versus the H/3 in- 
tensity are shown. Since the intensity of H/3 is proportional to the 
mass, the tendencies shown in these plots can be directly com- 
pared to the slopes of our scale-relations. Only the models that 
have an age, mass and metallicity for which our scale-relations 
predict EW >1 are plotted. Models with ages lower and greater 
than 4.5 Ma are plotted separately, because around that age many 
lines show a sudden change in their scale-relations slopes. Lines 
with slope equal to 1 and representative slopes from the scale- 
relations are shown as a reference for each emission line and 
age interval. The zero-point of these reference lines has been ad- 
justed for a better comparison with models. For [On] A 3727 
A and ages lower than 4.5 Ma the models follow a tendency 
with slope smaller than 1 , while for older ages the slope is near 1 . 
These tendencies coincide with the ones obtained with our scale- 
relations (see fig|5]l. The same occurs for [O ra] A 5007 A and 
[Nn] A 6584 A (Figs. [15] and [16]!. Moreover, the scatter in the 
plots, which is due to the different ages, metallicities and physi- 
cal parameters of the models, corresponds to the spread of values 



of the a coefficient of the scale-relations. These similar tenden- 
cies indicate that our scale-relations are robust against changes 
of density, inner radius and filling factor. 

5. Conclusions 

In this work we have studied the variation of the Hn regions 
spectrum with the ionizing cluster mass. For this task we have 
done photoionization models for several metallicities and ages 
in a cluster mass range from 10 3 to 10 7 M . From these models 
we obtain power-law scale-relations between the emission line 
intensities and the ionizing cluster mass. 

These power-law scale-relations are very useful tools since 
they allow to estimate in a simple way mass dependences of EW 
and line intensities ratios. Using such scale-relations we have 
checked that EW(Hfi) is practically independent of mass, a result 
which supports its use as an age indicator. On the other hand, 
EW^Oiii] 5007) has a mass dependence that, in spite of being 
small, implies that some caution is needed when using it as an 
age indicator. The mass dependence obtained for [O n] A 3727 
A implies that it is a good SFR indicator if it is calibrated with 
H ii regions older than 4.5 Ma, and used on objects dominated by 
such H ii regions. As the mass dependence of [O n] A 3727 A is 
similar to the mass dependence of [Nn] A 6584 A, the intensity 
ratio of these lines has a small mass dependence. This supports 
the use of that ratio as an abundance indicator. We also check 
that the R23 parameter is almost independent of mass, while S23 
and 77' have a non-trivial mass dependence. 

The estimation methods that use mass independent indica- 
tors can be calibrated with no consideration for the mass of the 
objects. If the indicator is mass dependent, however, the estima- 
tion method should be recalibrated taking into account the ioniz- 
ing cluster mass. For those methods, semi-empirical calibrations 
are, at least for now, more appropriate. 

An additional result obtained from our models is an intrinsic 
limit of about AE(B - V) = 0.1 in the accuracy of the extinction 
estimations carried out under the hypothesis assumed in synthe- 
sis models (case B and T e = 10 4 K). This error is diminished if 
tailored models are used. 

Finally, we have shown that the relevant collision lines have a 
non linear dependence with the cluster mass, and that, in general, 
not all the collisional lines have the same dependence (slopes in 
the log M - log L differs from line to line). A direct implication 
of this result is that different collisional lines in the emission 
line spectrum of a non-active galaxy have different contributions 
from different Hn regions and, hence, it depends on the Hn pop- 
ulation in the galaxy (and the underlying initial mass cluster dis- 
tribution) and that low mass clusters (< 10 4 M ) can contribute 
more than 50% of the luminosity of the line. In this situation, 
any attempt to infer physical properties from the emission line 
spectrum of a non-active galaxy by comparison with single neb- 
ulae photoionization models and without taking into account the 
IMF sampling effects may be a risky business. 
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Fig. 16. [Nil] A 6584 A / luminosity with respect to the Wfi luminosity for Stasihs ka et alj d2001l) models with ages lower (left) and 
higher (right) than 4.5 Ma. A solid line with slope equal to 1 and dotted lines with representative slopes are shown as references. 



